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Abstract 



High density phase transitions in a 4 dimensional Nambu-Jona-Lasinio model containing a single 
symmetry breaking order parameter coming from the fermion-antifermion condensates are researched 
' and expounded by means of both the gap equation and the effective potential approach. The phase 

^ , transitions are proven to be second order at a high temperature T; however at T = 0, they are first- 

QQ ' or second- order, depending on whether A/m(0), the ratio of the momentum cutoff A in the fermion 

fT^ ' loop integrals to the dynamical fermion mass m(0) at zero temperature, is less than 3.387 or not. 

The former condition can not be satisfied in some models. The discussions further show complete 
OO ' effectiveness of the critical analysis based on the gap equation for second order phase transitions 

' including determination of the condition of their occurrence. 

m ' 

O ■ PACS: ll.lO.Wx; ll.lO.Lm; ll.15.Pg; ll.30.Er 

' Keywords: Nambu-Jona-Lasinio model, symmetry restoring at high density, gap equation, effective po- 

tential, second and first order phase transition 



1 Introduction 

Nambu-Jona-Lasinio (NJL) model ^ is a good laboratory to research symmetry restoring phase transition 
at finite temperature T and finite chemical potential ^ [21. It could simulate the phase transitions in 
Quantum Chromo dynamics (QCD) [3] and is directly related to some dynamical schemes of electroweak 
symmetry breaking 01 |S] . It is known that in the NJL model with a single order parameter which comes 
from fermion-antifermion (//) condensates, the phase transitions at a high critical temperature Tc for 
fixed fi are second order <B . On the other hand, with regard to the phase transitions in this class of models 
at a high chemical potential for a fixed T, especially for T — + case, less deep-going researches were made 
and the conclusions failed to be made so simple and clear. Therefore, a more careful investigation of 
this problem is certainly interesting. On the high density phase transitions at T = in this class of 
NJL models, the authors of Ref. [7] derived the effective potential from the imaginary-time formalism 
of thermal field theory and briefly outlined the results. However, it seems that some more detailed 
demonstrations are necessary for understanding these results and in addition, some of the results there 
need to be clarified further. Actuated by the above motivates, in this paper, we will make a thorough 
examination of the phase transitions at high /i in this class of AD NJL model. The discussions will 
be made in the real-time formalism of thermal field theory and involve both T ^ and T — case. 
For second order phase transitions, we will make critical analysis of the order parameter of symmetry 
breaking based on both the gap equation and the zero temperature effective potential. For first order 
phase transitions we will explicitly give the equations determining the critical curves. The condition 
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which distinguishes between second and first order phase transitions will be indicated definitely. 
The Lagrangian of the model will be expressed by 

N N 

c{x) ^Y.^\x)irdM^) + lY,[^\x)M^)]\ (1) 

fc=l k=l 

where ipk{x) are the fermion fields with N "color" components and g is the four-fermion coupling constant. 
It is indicated that the four-fermion interactions in Eq.(l) could lead to only the fcrmion-antifermion con- 
densates {'ipip). The conclusions reached from Eq.(l) are essentially also applicable to more complicated 
models with a single order parameter coming from the condensates (iptp). The discussions will be made 
in the fermion bubble diagram approximation which is equivalent to the leading order of the ex- 
pansion. Since when D — A the Lagrangian (1) is not renormalizable, we will regard it only as a low 
energy effective field theory and examine its physical results. We will analyze second order high fi phase 
transitions based on the gap equation in Sect El examine second and first order high /i phase transitions 
by means of the zero temperature effective potential respectively In Sects. |3and01and in SectElcome to 
our conclusions. 



2 Second order high ^ phase transitions — gap equation analysis 

Assume that dynamical fermion mass could be generated by the fermion condensates induced by the 
four-fermion interactions in Eq.(l) at T = /i = and at some finite T and /i, then the corresponding 



mass term will spontaneously break the discrete chiral symmetry xd- i^kix) 75'0fc(a^) and the special 

v 

parities Vj: ipk{t, ■ ■ ■ ■ ■) — ^ j^ipk^t, • • • , —x^, ■ ■ = 1, 2, 3). After the gap equation at T = /i 
is substituted, the gap equation at T 7^ can be transformed to [S] 







-m2(0)ln 



A2 



m2(0) 



A2 

= -m^ In ( ^ + 1 



F3{T,fi,m) 



(2) 



where A is the four dimensional Euclidean momentum cutoff in the fermion loop integrals, m(0) and 
m = m(T, /i) are respectively the dynamical fermion mass at T = /i = and at T 7^ and/or n =/= and 



F3 (T, /i, m) = [h {y,-r) + l3{y, r)] 



with 



/3(.,T.) = i| 



dx 



\/ x"^ + exp(-\/a;^ + =F r) -f 1 



(3) 



(4) 



and y = m/T and r — ^/T . It is easy to verify that F^{T^ /i, m) is an increasing function of T and /i but 
a decreasing function of m, and from Eq.(2), we can assume that at some critical temperature Tc and 
critical chemical potential /ic, the dynamical fermion mass m becomes zero, thus the symmetries xd and 
Vj (j = 1, 2, 3) which are spontaneously broken at lower T and /i will be restored. Setting m = 0(jj — 0) 
in Eq.(2) we will obtain the critical equation 



-m2(0) In 



A2 



m2(0) 



1 



dx 



1 



(5) 



When taking the limit m — )■ and obtaining Eq.(5), we did not encounter any singularity, this means 
that TO may go to zero continuously thus Eq.(5) will be the critical equation of a second order phase 
transition. In fact, for a given high temperature T, substituting Eq.(5) with Tc replaced by T into Eq.(2), 
then using the high temperature expansion of F^{T, iJL,m)^, we can obtain the critical behavior of rn^ 
near /ic at high T 



= {A 
with fc = [i-cjT and 



^2^ cosh(prc 



■I 



In- 



A 
TV 



1 



+ 7-/i(T,m) 



when [i ^ and T is high, (6) 



Mr,A.) = 7C(3)(-f 



31C(5)(^)^ + 127C(7)(^)^ 



(7) 
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where C(s) (s — 3, 5, 7) are the Riemann zeta functions and 7 is the Euler constant. Equation (6) indicates 
that for high T, the phase transition at /Jc is second order. For proving the phase transitions at fj,c for 
low T, we first consider the limit case of T — > 0. In this case (r — > 00) , the critical equation (5) with 
He ^ ^^ and Tc T is reduced to 



= -m^(0)ln 



A' 



m2(0) 



/x^Q, when T ^ 0. 



Through changing the integral variable by z = {x^ + y^)^^^ we may rewrite 



F^lT, /i, m) = 2™^ 
with a — ii/m and obtain 

F3(T = O,/i,m) = 0(^-m) 



dz 



7 ^ 7 + i^Ol — > a) 



(8) 



(9) 



/x\/ /i^ — Tin? — rn? In 



(10) 



By using Eqs.(8) and (10), we may reduce the T limit of the gap equation (2) to 

m — m(0), when /i < m(0) 

when /X > to(0). 



2 2 

1^ ^ (m + \/Jj^~n?f 



(11) 



From Eq. (11) it can be proven that 



lim 



dr 



1 — 

Ofl 



In 



A 



= — 2-\/ /i^ — m? /m 
< O, when m(0) < fJ. < /ico 



2 , 2 

TO 



A^ 



— 7 

TO 



(/i + - m2)2 A 

2 ^2 



2 ™2 TO 



(12) 



and 



lim 



< 0, when to(0) < /i < /ico- 



(13) 



Eq.(12) shows that in the T ^ limit, when = to(0), dm/dfi = and when /i = /ico where m = 0, 
dm/dfi = —00. From these results and Eq.(13), we can deduce that the m — ^ curve at T = is concave 
downward in the region to(0) < fi < HcO- Near the critical point ~ /Xco, we find out from Eq.(ll) that 



TO^ ~ {hIq - [1^)1 [ln(A/2/x) - 1/2] , when T ^ 0. 



(14) 



Eq.(14) indicates that in the D — \ NJL model with a single order parameter, when T = 0, the same as 
when T is high, the symmetry restoration phase transition at [i^ could be second order. This is different 
from the high density phase transitions when T = in Z) = 2 and D = 3 GN model where they are only 
first order HJinilini. 

However, for consistence of Eq.(14) i.e. when /i < /icOj > 0, we must have ln(A/2/i) — 1/2 ~ 
ln(A/2/ico) — 1/2 > and this implies that 



Mco < AV4e 

which is the condition in which a second order phase transition could occur at T = 0. 



(15) 



3 second order phase transitions — effective potential analysis 

The above conclusion that when the condition (15) is satisfied the high density phase transition at T = 
is second order reached only by the gap equation analysis can be verified by an effective potential analysis. 
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It has been proven JT] that the extreme value condition of the thermal effective potential V^pj:{T, ii,m) 



IS 



with 



dv}^l.{T,fi,m) _ m 



dm 



27r2 



■In 



dm 



1 



m2(0) 



In 



A2 



m2(0) 



1 



F3{T,^i,m) 



(16) 



(17) 



Hence Eq.(16) just corresponds to the gap equation (2) multiplied by m. From Eqs. (17) and (10) we can 
find out the effective potential V^^^j.-{T — 0,ii,m) at T = satisfying the condition V"j'^j (T — 0, /i, m 



0) = expressed by 



K';}(r = 0,M,m) 



1 
2P 



A* 
"8" 



In 1 



m 
A^ 



m^(0)TO^ 



In 



A^ 



m^(0) 



+^ + 9{fi - m) 



^ [L^ — vn? — — m?) ^ 



^^In. 



(18) 

Eq.(18) is essentially identical to the effective potential in Ref. derived from the imaginary-time 
thermal field theory, except the extra "normalization" condition Vg^|(T = 0, ^,to = 0) = 0. From Eqs. 

(10) and (17) we can directly obtain the extreme value condition of yj^|(r = 0,/i,m) 



ayi;)(T = o,/i,m) 

dm 



In 



A^ 

m 



m^(0) 



In 



A^ 



m^(0) 



li'- — m'- 



'■ In ' 



/.J? — m^ 



= 



6'(/x — to) 



(19) 



and 



5Vj,'|(r = 0,/x,m) 



dm^ 



1 

2^ 



fTO^ln 



A^ 
m 



A" 



TO^(O) 



In 



A^ 
^^7(0) 



+6'(Ai - to) 



jj? — — Sto^ In ■ 



+ -yZ/i^ — TO^ 



(20) 



(4') 

Variation of V;y}(r 0, /i, to) as /i increases may be discussed by Eqs. (18)-(20). 

First we note that when /i = 0, for any K^/m'^{Q) ^ 0, Vg'^|(r = G,/i = 0, m) will have a maximum 
point m = and a minimum point m = to(G) and this shows spontaneous symmetry breaking at 
T = /i = 0. 

If the symmetries will be restored at high /i through second order phase transitions, then the extreme 
value equation of v}^j{T = 0, /x, to) with m ^ 0, i.e. the T — > limit of the gap equation (2), should have 

the only solution and it must correspond to a minimum point of ^^^^(T = 0, fi, m). We will examine this 
case in the condition expressed by Eq.(15). This condition, by Eq.(8), is 



m2(0) 



In 



A2 



to2(0) 



< 



A2 

4e 



which amounts to A/to(0) > 3.387 and will lead to ^co > m{0). Variation of Fg^|(r — 0,/i, m) as ^ 
increases is as follows. 

1) < /i < to(0). In this case, V^^j{T ~ G,//,to) will have a maximum point to = 0(when m < fi) 
and a minimum point to = TO(0)(whcn to > /^). It can be proven that the other extreme value equation 
contained in Eq. (19) when m < ji 



iTO^(0)ln 



A2 



m2(0) 



■In 



(^ + \J — to2)2 



+ /i2 — to2 



(21) 
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has no solution for m < ^ < m(0) when i In 
then Eq.(21) can be changed into 



TO^(O) 



+ 1 



> 1. In fact, if we set fj, = am(0) and m = /3/x, 



A2 



m2(0) 



(3^a\ AVm2(0) + a2/32 
— - — In ■ 



a2(l + ^1^)2 



(22) 



It is easy to check by numerical solution that Eq.(22) has no solution for a < 1 and /3 < 1 when 



iln 



(0) 



1 



> 1 or A2/m2(0) > - 1. Our prerequisite A/m(0) > 3.387 or A^/m'^{0) > 15.11 
apparently satisfies this condition. Thus Eq.(21) does have no solution. As a result, m = m(0) becomes 
the only minimum point of V^jj{T = 0, /x, m) and we will have the same spontaneous symmetry breaking 
as one at T = = 0. 

2) m(0) < fJ. < 1-J.co- In this case, when m > ji, y,'jj(T = 0,fi,m) has no extreme value point and 

when m < jj,, V^jj:{T = 0,/i,m) will have a maximum point m = and a minimum point nii which 
is determined by Eq. (21). Noting that in present case Eq.(21) may have the solution m = m\ when 



:ln 



to2(0) 



+ 1 



> 1 and 



dw? 



A2 



2 A2 + ; 



> 0, when < ^,co. 



In fact, Eq.(21) has the solutions 



and 



mi = m(0),when fi = m{0) 



mi < m(0),when fj, > m(0). 



(23) 



(24) 



noting that equation (21) obeyed by mi is just the second formula in Eq.(ll), hence we can obtain from 
Eq.(12) that dmi/dfi < 0. Equation (23) indicates that when fi = m(0), we will go back to the same case 
as n < m(0) and Eqs. (21) and (24) show that as fi increases from m(0) further, the minimum point toi 
will become smaller and smaller and finally it will go to zero continuously at a critical chemical potential. 
Correspondingly, the global minimum of V^^^j.-{T = 0,/x, m) at m = mi will go up as fj, increases from a 
negative value to zero, since it can be proven that 



<)(T: 



0,/U,mi) 



127r2 



3toi 



+ 4/x3-(/x2-m?)3/2-3M' 



> 0. 



The critical chemical potential can be determined by taking m = in Eq. (21) and the result is precisely 
the /XcO given by Eq.(8). Since Eq.(21) is exactly the second formula in Eq.(ll) and the critical behavior 
(14) of m2 follows. Equation (8) expresses a critical curve C2 of second order phase transitions in /x — m(0) 
plane. 

3) /X = jjLco- Now V^^^^{T = 0, /Lt, m) has the only extreme value point m = and it is found out that 



S"yJ/}(r = 0,M,m) 



dm" 



m=0 



3 
2^ 

3 



In 



0, 
A2 



- 1 



^0 



when n = 2, 3, 5 
when n = 4 

when n = 6 



(25) 



Equation (25) implies that in the condition \n{A'^ /Afilo) — 1 > or ^"^q < A2/4e, m = will be the only 
minimum point of V"Jjj(T = 0, /z, m) and the broken symmetries will be restored. 

4) /X > /Xco- In present case, it can be proven that the extreme value equation (21) has no solution when 

m < ji. This can be seen most simply from Eq.(14) which is the approximation of Eq.(21) when /i ~ /Xco 
and m ~ 0. More rigorously, we can set /i = 7/(Xco, m = /3/x = PjjJcO, then Eq.(21) can be changed into 



In 



A2 



21n[7(l + v/r^)] 



/J2 



1 



(26) 
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It may be checked that Eq.(26) has no solution with 7 > 1 (/i > /7,co) and /3 < 1 (m < /i) when 
Is^Q < A^/4e. Consequently, ^^^^(T = 0,yLt,m) will have the only minimum point to = in this case and 
this fact further indicates restoration of the symmetries which were broken at T = and jj, < /XcO through 
a second order phase transition when Eq.(15) is satisfied. 

Up to now we have proven that the phase transitions at jJc are second order for a high T and for 
T = when /z^q < A^/4e. Considering these results and the fact that Eq.(5) is a critical equation of 
second order phase transitions for any finite Tc and ^c, we can conclude that the phase transitions at fXc 
for any T including low T are second order when fM^Q < A^/4e. 



4 First order phase transition at T = 

A first order phase transition could generally occur in the case when the extreme value equation of 
V^'j|(T = 0, jjL, to) has two or more solutions with to ^ 0. We will prove that is the case when /x^g > A^/4e 
and fi > jico ■ First we indicate that when /i > /ico , to = is always a minimum point ofV^jjiT = 0,/i,m) 
and when /x^q > A^/4e , Eq.(26) will have the solutions with 7 > 1 and /3 < 1, i.e. the extreme value 
equation (21) will always have solutions for jj, > jjLco > A/2e^/^ and m < jjL. This conclusion can also be 
checked by examining Eq.(14). Since in this case to = is a minimum point, the solutions of Eq.(21) 
must correspond to at least a maximum and a minimum point of V^^^j:(T = 0, /i, to) and this could lead 
to a first order phase transition. Assume the maximum point is TO2 and the minimum point is TOi, then 
a first order phase transition curve should be determined by the equations 



0,n,m = mi) = v}fl{T: 
0,IJ.,m) 



0,li,m = 0)= 0, 



dm 



m=mi7^0 



and 



a2l/W(T = 0,M,m) 



dm? 



> 



(27) 
(28) 

(29) 



m=0 



whose explicit or equivalent forms are 



mf In ( — 2 + 1 



TOt 



■A2 



A4ln(l + ^ ) - 4to?Mc+ oM' + ^(M-mi) 



2/iTOiY^/Lt2 - ml - ^nilJ? - TOi)^/^ - TOi In 



TOT 



= 0, 



!^ln(| + l)+.(,-TO.) 



TOi 



(30) 



(31) 



where toi ^ and fj, > Hco- When n < mi, from Eqs.(31) and (8) we obtain toi = to(0). It may be pointed 

out that the condition /i^Q > A^/4e means that A/m(0) < 3.387 thus we can have cither /ico > to,(0) 
or /Xco < m{0), so even if > /XcO) it is still possible for /x < to(0). Substituting /x < toi = to(0) into 
Eq.(30), we are led to that 



3^^ 



to4(0) 



ln(a + 1) - a + In ( 1 + - 



A2 



to2(0)' 



(32) 



which expresses a first order phase transition curve Ci in the fj, — to(0) plane. Eq.(32) is valid only if 
H < nil = TO,(0), this will lead to the constraint ln(a + 1) - a + In (l + ^) < 4/3 thus A/to(0) < 2.21. 
So the curve Ci must be in the region of to(0) > A/2.21. 
When /i > to(0), Eqs. (30)-(31) can be reduced to 



- - mlf/' = 2to?m^o - ^'rnj - AHn (l + ^) , 



(33) 
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7 + 



2 

'1- 



(34) 



— TO- 



2 (^+ ^/i2 - mj) 



Equations (33) and (34) together with toi ^ and /i > /i^o will be the equations to determine the first 
order phase transition curve C[. When /i to(0), we may obtain from Eq.(34) that nii to(0), as a 
result, equation (33) becomes 



which coincides with Eq.(32) of the curve Ci at the point /i = to(0). This means that the two first order 
phase transition curves C[ and Ci meet at the point — to(0). On the other hand, we note that fi —^ /ico 
from /J, > /ico and toi — > is a limiting solution of Eqs.(33) and (34) and that the curve C[ should be in 
the region //^q > A'^/4e, so it must intersect the second order phase transition curve C2 at //^q — h? /Ae. 
Consequently, in the /j, — to(0) plane the curve C( will start from the point /j,^q = A2/4e, extend itself in 
the region of /i > /^co and end at /i = to(0) where it meets the curve Ci. The curve C'l must be in the 
region of /i > /Xco, since for the opposite case of /i < /ico, it may be seen from Eqs. (20)and (25) that 
when /i^Q > A2/4e, m = will be a maximum point of V^^^^{T = 0,^,to), and this only corresponds to 
the phase of symmetries being broken. The above discussions show that as far as symmetry restoring 
phase transition is concerned, fi^ — /i^Q — /Ae will be a tricritical point in the /i — to(0) plane. We have 
noted that the similar results were obtained in Ref. 7 , however, no detained demonstration leading to 
the results, especially no distinction and relation between the first order phase transition curves Ci and 
C'l were given there. 



5 Conclusions 

In this paper, we have researched the phase transitions at high density in a 4D NJL model with a single 
order parameter of symmetry breaking coming from fermion-antifermion condensates by means of the 
gap equation and the zero temperature effective potential. We have proven that for high T the symmetry 
restoring phase transitions are always second order; and for T = 0, depending on whether /x^g is bigger 
than h? /Ae or not, or equivalently, whether A/m(0) is less than 3.387 or not, the phase transitions will 
be first- or second- order. Hence //^q — K'^ /Ae is a tricritical point of high density phase transitions at 
T = 0. We have also further deduced that the phase transitions at fjLc for any T are second order when 
A/to(0) > 3.387. 

We note that in this class of models, first order phase transitions occur only if the dynamical fermion 
mass to(0) at T = is close to and has the same order of magnitude as the momentum cutoff A of 
the loop integrals. However, this condition is generally not natural or can not be satisfied in some low 
energy effective theories of NJL-form. For instance, in the top-quark condensate scheme of electroweak 
symmetry breaking 21 j A/to(0) must be so large that up to ~ 10^^ owing to the constraint from the 
standard electroweak model. Even if in the four-generation fermion extension of such scheme 5 where 
A/to(0) could go down greatly, it must still be bigger than 10. Therefore, in some physical 4D NJL 
models, it can be assumed that first order phase transition in fact does not occur and even if at T = 0, 
there are only second order phase transitions. 

It is seen from the above discussions that as far as analysis of a second order phase transition is 
concerned, the gap equation approach and the effective potential approach have the same effectiveness, 
including determination of the occurrence condition of a second order phase transition. However, the 
critical analysis of the order parameter m based on the gap equation is apparently more simple and direct 
than the ones based on the effective potential. 

The total conclusions in this paper are obtained only in the model (1) with a single order parameter 
coming from the fermion-antifermion condensates (^/'?/') and they could have substantial change when 
one transfers to the model with both the fermion-antifermion condensates ('0V') ^^nd the fermion-fermion 
condensates (V'V')- The latter model may be a better simulation to QCD and deserves to be researched 
further. 

The author is grateful to Professor K.G. Klimenko for a private communication which brought the 
author's attention to the relevant work in Ref.(2j- 
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